THE METHOD OF BRACKETS. PART 2: EXAMPLES 
AND APPLICATIONS 



IVAN GONZALEZ, VICTOR H. MOLL, AND ARMIN STRAUB 

Abstract. A new heuristic method for the evaluation of definite 
integrals is presented. This method of brackets has its origin in 
methods developed for the evaluation of Feynman diagrams. The 
operational rules are described and the method is illustrated with 
several examples. The method of brackets reduces the evaluation 
of a large class of definite integrals to the solution of a linear system 
of equations. 



1. Introduction 

The method of brackets presented here provides a method for the eval- 
uation of a large class of definite integrals. The ideas were originally 
presented in [6] in the context of integrals arising from Feynman dia- 
grams. A complete description of the operational rules of the method 
together with a variety of examples was first discussed in [5]. 

The method is quite simple to work with and many of the entries 
from the classical table of integrals [7\ can be derived using this method. 
The basic idea is to introduce the formal symbol (a), called a bracket, 
which represents the divergent integral 



/ x a ~ l dx 



The formal rules for operating with these brackets are described in 
Section [2] and their justification (especially of the heuristic Rule 12. 3p is 
work-in-progress. In particular, convergence issues are ignored at the 
moment. Roughly, each integral generates a linear system of equations 
and for each choice of free variables the method yields a series with the 
free variables as summation indices. A heuristic rule states that those 
converging in a common region give the desired evaluation. 
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Section [3] illustrates the method by evaluating the Laplace trans- 
form of the Bessel function J v (x). In this example, the two resulting 
series converge in different regions and are analytic continuations of 
each other. This is a general phenomenon which is used in Section 
[5] to produce an explicit analytic continuation of the hypergeometric 
function q+ iF q (x). Section @] presents the evaluation of a family of in- 
tegrals C n appearing in Statistical Mechanics. These were introduced 
in jl] as a toy model and their physical interpretation was discovered 
later. The method of brackets is employed here to evaluate the first 
four values, the only known cases (an expression for the next value C5 
in terms of a double hypergeometric series is possible but is not given 
here). The last section employs the method of brackets to resolve a 
Feynman diagram. 



2. The method of brackets 

The method of brackets discussed in this paper is based on the as- 
signment of the formal symbol (a) to the divergent integral (11. ip . 

Example 2.1. If / is given by the formal power series 

00 

n=0 

then the improper integral of / over the positive real axis is formally 
written as the bracket series 

POO 

(2.1) / f(x)dx = J2a n (an + P). 

Jo 



Here, and in the sequel, Ylm ls use d as a shorthand for 



n=0" 



Formal rules for operating with brackets are described next. In par- 
ticular, Rule 12.21 describes how to evaluate a bracket series such as the 
one appearing in (12.11) . To this end, it is useful to introduce the symbol 

(-l) n 

(2.2) 0, - 



r " r(n + i)' 

which is called the indicator of n. 

Example 2.2. The gamma function has the bracket expansion 

/■oo 

(2.3) r(a) = / x a - 1 e- x dx = y2 ( l ) n(n + a). 

Jo 



METHOD OF BRACKETS 



3 



Rule 2.1. The bracket expansion 
(2.4) 

1 j. mi „m. (a + m 1 + --- + m r ) 



(ai + a 2 H h a r )° ' r(a) 

v ' mi,..., rrir v ' 

holds. Here 4> mi ,...,m r is a shorthand notation for the product <p mi ■ ■ ■ 4> mr . 
If there is no possibility of confusion this will be further abridged as 
4>{m}- The notation J2{ m } * s ^° ^ e understood likewise. 

Rule 2.2. A series of brackets is assigned a value according to 

(2.5) V <f> n f{n) (an + b) = -L/(„*)r(-n*), 

a 

n 1 1 

where n* is the solution of the equation an + 6 = 0. Observe that this 
might result in the replacing of the index n, initially a nonnegative 
integer, by a complex number n*. 

Similarly, a higher dimensional bracket series, that is, 

yj 0{n}/(^i, • • • , n r ) (anUi + • • • a lr n r + c\) ■ ■ ■ (a r \n\ + • • • a rr n r + c r ) 

{n} 

is assigned the value 

(2-6) ]d^(Aj\ fin ^ ' ' ' ' n *) r ^ n i) • • • r ^<)' 

where A is the matrix of coefficients (a^) and (n*) is the solution of 
the linear system obtained by the vanishing of the brackets. The value 
is not defined if the matrix A is not invertible. 

Rule 2.3. In the case where a higher dimensional series has more sum- 
mation indices than brackets, the appropriate number of free variables 
is chosen among the indices. For each such choice, Rule 12.21 yields a 
series. Those converging in a common region are added to evaluate the 
desired integral. 

3. An example from Gradshteyn and Ryzhik 

The second author is involved in a long term project of providing 
proofs of all the entries from the classical table of integrals by Grad- 
shteyn and Ryzhik |7j. The proofs can be found at: 



fli • • • Qj 



http : // www . math . tulane . edu/~ vhm/Table . html 



In this section the method of brackets is illustrated to find 



(3.1) / x u e~ ax J v ((3x) dx 



v -ax j ( n\ i _ V r ) \ 2' 



^(a 2 + (3 2 Y +1 / 2 
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which is entry 6.623.1 of [7]. Here 

-l) k (x/2) 2k+v 



(3.2) J v {x) = 



, o k\ T(k + v + 1) 

is the Bessel function of order v. To this end, the integrand is expanded 
as 



(3.3) e-°*J u (/3x) = ^ £ MaxTj ^ V{k + u + l) 

" \2> „n+2k+2v 

T(k + U+1) 



= ^2<Pk.n— " 7T J ' 

k,n 

so as to obtain the bracket series 

poo qJI/ P_\2k+v 

(3.4) / e- ax J u (f3x)dx = V <j> k>n — ^ — (n + 2k + 2v + 1) . 

The evaluation of this double sum by the method of brackets produces 
two series corresponding to using either k or n as the free variable when 
applying Rule 12.21 

The index k is free. Choosing k as the free variable when applying 
Rule 12. 21 to (13. 4p . yields n* = —2k — 2v — \ and thus the resulting series 

a -2k-2v-lf§_\2k+v 

^ r(fc + z/ + i) 



.^^ £(2^+1) (v+\ 

a ™ r>TT) lFo 



a 2 



The right-hand side employs the usual notation for the hypergeometric 
function 



(3.6) p F q 



oo 

d\ , ... , dp 



h,...,b, 



X 



\- («i)n • • • (a P ) w x n 



where (a) n = r ^^ n - > is the Pochhammer symbol. Note that the \Fq in 
(13.51) converges provided \f3\ < \a\. In this case, the standard identity 
iFq(o\x) = (1 — x)~ a together with the duplication formula for the V 
function shows that the series in (13. 5p is indeed equal to the right-hand 
side of (|3II1). 
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The index n is free. In this second case, the linear system in Rule 
has determinant 2 and yields k* = —n/2 — v — 1/2. This gives 



1 ^ a r 



(3.7) ly^fh^^ r(n/2 + i/ + l/2). 

v ) T -n/2 + 1/2 v ' ' 1 

This series now converges provided that > | a | in which case it again 
sums to the right-hand side of ( 13. II) . 

Note. T/iis z's t/ie typical behavior of the method of brackets. The dif- 
ferent choices of indices as free variables give representations of the 
solution valid in different regions. Each of these is an analytic contin- 
uation of the other ones. 



4. Integrals of the Ising class 
In this section the method of brackets is used to discuss the integral 
4 f°° f°° 1 d Ul du n 



(4.1) C n 



nl Jo Jo (sr n („ . U, \\ 2 u i u r 



(e; =1 k + 1/S))' 



This family was introduced in [4] as a caricature of the Ising suscepti- 
bility integrals 
(4.2) 

D _ ± f°° f°° -Q fui - u A 2 1 r///, da 



111 •/(] ./(i " V", - ".// (sr !i 1,, i / w '"J "» 



Actually, the integrals C n appear naturally in the analysis of certain 
amplitude transforms [10] . The first few values are given by 

(4.3) C x = 2, C 2 = 1, C 3 = L_ 3 (2), C 4 = ^C(3). 

Here, Lp is the Dirichlet L-function. In this case, 

("> ^ = E ((STTT? " (STT^J ' 

No analytic expression for C n is known for n > 5. Similarly, 
(4.5) 

1 47T 2 47T 2 1 7 

£>! = 2, D 2 = -, D 3 = 8 + — - 27L_ 3 (2), D 4 = — - - - —((3) 
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are given in p£j- High precision numerical evaluation and PSLQ exper- 
iments have further produced the conjecture 

(4.6) 

D 5 = 42 - 1984Li 4 (|) + ^vr 4 - 74C(3) - 1272C(3) In 2 + 40vr 2 ln 2 2 

- — 7T 3 + y 7T 2 In 2 + 88 In 4 2 + 464 In 2 2 - 40 In 2. 
The integral C n is the special case k = 1 of the family 

4 r r i d Ul d Ur 



that also gives the moments of powers of the Bessel function K via 

on— fc+1 on— fc+l Z" 00 

The values 

(4.9) c^W^Y c,,,-^^^) 



2 ; ' AK 4 r (| + 1) ' 

as well as the recursion 

(4.10) (k + l) 4 c 3 , fc - 2(5k 2 + 20k + 21)c 3 , fc+2 + 9c 3 , fc+4 = 

with initial data 
(4.11) 

3« 3 a 4vr 5 2 

C3,0 — 7^— , C 3) i — ~L,- 3 {;&), C 3 2 — — ——, C 3 3 — iv_ 3 ^J - -, 

327r 4 967T 9a 3 

where a = 2~ 2 / 3 r 6 ( 1 r) are given in pQ and [3]. 

The evaluation of these integrals presented in the literature usually 
begins with the introduction of spherical coordinates. This reduces the 
dimension of C n by two and immediately gives the values of C\ and 
C*2. The evaluation of C 3 is reduced to the logarithmic integral 

2 f°° ln(l + x) dx 



( 4 - 12 ) Cs o , 

O Jq X — r 3? "T" -L 

Its value is obtained by the change of variables x — > | — 1 followed by 
an expansion of the integrand. A systematic discussion of these type 
of logarithmic integrals is provided in [9]. The value of C4 is obtained 
via the double integral representation 

1 f°° f°° \n(l + x + y) dxdy 



1 r°° r c 

(413) c <=iL L 



[1 + x + y)(l + 1/x + l/y) - 1 x y 
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Moreover, the limiting behavior 

(4.14) lim C n = 2e" 27 

n—too 

was established in jl]. 

In this section the method of brackets is used to obtain the expres- 
sions for C 2 , C 3 , and C 4 described above. An advantage of this method 
is that it systematically gives an analytic expression for these integrals. 
When applied to C 5 , the method produces a double series representa- 
tion which is not discussed here. 

4.1. Evaluation of C^fc. The numbers C2.k are given by 

dx dy 



POO POO 

(4.15) C 2 , k = 2 / 

Jo Jo 



lo Jo xy (x + 1/x + y + l/y) k+1 
A direct application of the method of brackets, by applying Rule 12.11 
to the integrand as in (14.151) . results in a bracket expansion involving 
a 4-fold sum and 3 brackets. Rules 12.21 and 12.31 translates this into 
a collection of series with 4 — 3 = 1 summation indices. However, it 
is generally desirable to minimize the final number of summations by 
reducing the number of sums and increasing the number of brackets. 
In this example this is achieved by writing 

(xy) dx dy 



C 



Jo (% 2 y + y + xy 2 + x) k+l 

(xy) k dx dy 



2, A' 



2 

lo Jo (xy[x + y} + [x + y]) k+1 ' 

In the evaluation of these expressions, the term (x + y) must be ex- 
panded at the last step. The method of brackets now yields 

1 



_ = (k x m y m u + v) m + m (k + I + n, + n 2 ) 



( Xy [x + y ] + [x + y ]) k+1 ^ 2 r <* + ^ 

and the expansion of the term (x + y) gives 



713, "-4 



Replacing in the integral produces the bracket expansion 

n _ o a (k + l+ n 1 +n 2 ) (-n x - n 2 + n 3 + n 4 ) 

<^2,k — * /_,<P{n} 7777 ~T\ 777 s 

^— ' T(fc + 1) r(-ni — n 2 ) 

x (k + 1 + m + n 3 ) (k + 1 + rii + n 4 ) . 

The value of this formal sum is now obtained by solving the linear 
system k + 1 + rii + n 2 = 0, —n x — n 2 + n 3 + n 4 = 0, k + 1 + rii + n 3 = 0, 
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and k + l + ni+n4 = coming from the vanishing of brackets. This 
system has determinant 2 and its unique solution is n\ = n\ = n| = 
n * = _*±1. It follows that 

U]6) r T(-nl)T(-n*)T(-nl)T(-nl) _ T (^) 4 
1 ' ' 2 ' k T(k + l)r {-n\ - n* 2 ) T(k + l) 2 ' 

Note that, upon employing Legendre's duplication formula for the T 
function, this evaluation is equivalent to (14. 9p . In particular, this con- 
firms the value C 2 = C 2i i = 1 in (14. 3ft . 

Remark 4.1. The evaluation 

(4.17) 

x a ~ l yP~ 1 dx dy 



OO POO 



lo (x + l/x + y + l/y) k+1 



T(k + 1) 2 

that generalizes C 2j fc is obtained as a bonus. Similarly, 

(A 1^ t f R\ o f°° f°° ^V-'dxdy 
(4.18) J rs {a,p) = 2 



Jo (x + y) r (xy + l) s 

Y ^ -r+q+/3 j p ^ 2s+r-q~/3 \ p p ^ r-a+/3 ^ 

= W) ' 

Note that C 2>fc (a, /3) = J fc+ i, fc+ i(a + fc + 1, /3 + fc + 1). 

Remark 4.2. The Ising susceptibility integral _D 2 , see ( 14. 2fl . is obtained 
directly from the expression for J r s given above. Indeed, 



OO POO 



(4, 9) C! = 2 y o k (*-^ {x *££ +1 ? 

= 2 ( J 4 , 2 (4, 2) - 2 J 4>2 (3, 3) + J 4 , 2 (2, 4)) 
_ 1 
~ 3' 

This agrees with ( 14. 5ft . This technique also yields the generalization 
(4.20) D 2M )^r ff^V --V- 1 ^ 2 

7o Jo v^ + yy (x + i/x + y + i/y) 2 

(6-a)(6 + a)(2 + (6 - a) 2 )vr 2 
12(cos(a7r) — cos(/97r)) 

with limiting case Z) 2 (a,a;) = ~ . T" ■> • 

z V ' / 3 sin(o7r) 
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4.2. Evaluation of Cz,h- Next, consider the integral 
(4.21) 

2 r°° r°° r°° dxdydz 

3 Jo J J X yz (x + 1/x + y + l/y + z + l/z) k+1 
2 ,oo ,oo p (xyz) k dxdydz 



3 7o ./o </o (xyz (x + y) + z (x + y) + xyz 2 + xy) k+1 

The second form of the integrand is motivated by the desire to to 
minimize the number of sums and to maximize the number of brackets 
in the expansion. The denominator is now expanded as 

V d> {ny (xyT Mn *z^+ n * +2n s (x + yP +n2 (fc + 1 + ni + n2 + n3 + n4) 
W T{k + 1) 

and further expanding (x + y) ni+ ™ 2 as 

(-ni -n 2 + n 5 + n e ) 



r(-ni - n 2 ; 



produces a complete bracket expansion of the integrand of C^. Inte- 
gration then yields 

(4.22) 

n _ 2 1 , (—ni -n 2 + n 5 + n 6 ) 

<^3,k — -TT 2>^{™} FT ^ 

{"} 

x (k + 1 + rix + n 2 + n 3 + n 4 ) (k + 1 + rix + n 3 + n 4 + n 5 ) 
x (A; + 1 + rii + n 3 + n 4 + n 6 ) (fc + 1 + n x + n 2 + 2n 3 ) . 

This expression is regularized by replacing the bracket (k + 1 + n\ + 
n 2 + 2n3) with (k + l + rii + n 2 + 2n 3 + e) with the intent of letting e — > 0. 
(This corresponds to multiplying the initial integrand with z e ; however, 
note that many other regularizations are possible and eventually lead to 
Theorem l4.3l It will become clear shortly, see f)4.24p . why regularizing is 
necessary.) The method of brackets now gives a set of series expansions 
obtained by the vanishing of the five brackets in (14.221) . The solution 
of the corresponding linear system (which has determinant 2) leaves 
one free index and produces the integral as a series in this variable. Of 
the six possible free indices, only n 3 and n± produce convergent series 
(more specifically, for each free index one obtains a hypergeometric 
series %F 2 times an expression free of the index; for the indices n 3 ,n4 
the argument of this 3 F 2 is ~ while otherwise it is 4.) The heuristic 
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Rule 12.31 states that their sum yields the value of the integral: 
(4-23) C 3 , fc = - lim - -^-(f^) + fk,n{-e)) 

n=0 

where 

< 4 - 24 > = S^nfrr^ • 

Observe that the terms fk, n ( e ) are contributed by the index n 3 while 
the terms fk,n{—z) come from the index n±. At e = 0, each of them has 
a simple pole. Consequently, the even combination fk, n ( e ) + fk,n(~ e ) 
has no pole at e = 0. Using the expansions 

(4.25) F(x + e) = r(a;)(l + if>(x)e) + 0(e 2 ), 
for i^O, —1, —2, . . ., as well as 

(4.26) T(-n + e) = f-+V(n + l)] + 0(e), 

n! \e / 

for n = 0, 1, 2, . . ., provides the next result. 
Theorem 4.3. The integrals C^k o,re given by 

71=0 7 

In particular, for = 1 
(4-27) C7 3 = -X: 7^ mn + 2) - V(n + 1)) . 

n=0 ^ 

The evaluation of this sum using Mathematica 7 yields a large col- 
lection of special values of (poly-) logarithms. After simplifications, it 
yields C 3 = L_ 3 (2) as in ( jjj) . 

Remark 4.4. An extension of Theorem 14.31 is presented next: 
(4.28) 

x a ~ l y"~ 1 z 1 ~ 1 dx dy dz 



C 3>k (a,(3,y) 



for 7 = 0, is given by 



OO POO POO 



(x + l/x + y + l/y + z + l/z) k+1 ' 



i r fn + k+1±a±/3 ) / i 

M ' 2 M^ ri + 1 )_i v ,( n+ fc±l±aM) + ^( 2n + jfe + l) 



fc! ^ N) 2 r(2n + fc + i) V 2 



where the notation r(n+ = r(n+ fc+1 t a+ ^ )I> + 



2 



as well as V(n + k+1±a± ? ) = ^(n + k + 1 + a +P ) + V(n + fc+1 + Q ~^ ) + • 
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is employed. Similar expressions can be given for other integral values 
of 7. In the case where 7 is not integral, C 3) fc(a, (3, 7) can be written 
as a sum of two 3 -F 2 's with T factors. The symmetry of C 3i fc(a, 0, 7) in 
a, (3, 7, shows that this can be done if at least one of these arguments 
is nonintegral. 

4.3. Evaluation of C 4 . The last example discussed here is 
1 f°° f°° f°° f 00 dxdydzdw 



2 ■ 



6 Jo Jo Jo Jo xyzw (x + 1/x + y + 1/y + z + 1/z + w + 1/w) 

To minimize the number of sums and to maximize the number of brack- 
ets this is rewritten as 

1 f°° r°° r°° r°° x 1+e y l+e z 1+e w l+e dxdydzdw 

6 Jo Jo Jo Jo [Axyzw(x + y) + zw(x + y) + xyzw(z + w) + xy{z + w)\ 

with the intent of letting e — > and A — > 1. As in the case of C 3j fc, 
the regulator parameter e is introduced to cure the divergence of the 
resulting expressions. Similarly, the parameter A is employed to divide 
the resulting sums into convergence groups according to the heuristic 
Rule 12.31 The denominator expands as 

w 

x (x + y) ni+n2 (z + w ) n3+ni (2 + m + n 2 + n 3 + n A ) . 

As before, 
and 

(* + = £ <-"»-'■'+ "'+-«) . 

' r(-n 3 - n 4 ) 

These expansions of the integrand yield the bracket series 
(4.29) 

- ^0{n} A ni (2 + m + n 2 + n 3 + ti 4 ) 



6 

{n} 



(—Til - n 2 + n 5 + rie) (— n 3 - n 4 + n 7 + n 8 ) 
r(-m-n 2 ) r(-n 3 -n 4 ) 
x (2 + e + rii + n 3 + n 4 + n 5 ) (2 + e + rii + n 3 + n 4 + n 6 ) 

x (2 + e + rii + n 2 + n 3 + n 7 ) (2 + e + n x + n 2 + n 3 + n 8 ) . 
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The evaluation of this bracket series by Rules 12.21 and 12.31 yields hy- 
pergeometric series with arguments A (ni, n 2 , n 5 , or n 6 chosen as the 
free index) and 1/A (n 3 , n 4 , nj, or n 8 chosen as the free index). Either 
combination produces an expression for the integral C4. Taking those 
with argument A (the indices and tiq yield the same series; however, 
it is only taken into account once) gives 

(4.30) 



-A-T 2 ( e )r 2 (l-e) 



A c 



+ 



1 + 2e 

A~ e 



1 -2e 

As e — > 0, the limiting value is 



2F1 
2F1 



l + e 

±-6,1 



— e 



.4 



A - 2 2 F 1 



1.1 

3 
2 



.4 



(4.31) 



1 

24 



In 2 A In 



1 + V 7 ^ 1 



.4 



+ 



3VA 
In A 
6y/A 



u 3 (Va) - u 3 (-Va) 
u 2 {Va)-u 2 {-Va) 



Finally, the value of C4 is obtained by taking A — > 1: 
(4.32) C 4 = 

This agrees with (14. 3p . 



i[Li 3 (l)-Li 3 (-l)] = ^C(3). 



5. Analytic continuation of hypergeometric functions 
The hypergeometric function p F q , defined by the series 



(5.1) 



F„(x) 



v ± 1 



1 dp 



ai, 
bx,...,b q 



x 



E 

ra=0 



(fll)n 



(h) n ■ ■ ■ (bq)n Tl\ 



converges for all x G C if p < q + 1 and for \x\ < 1 if p = q + 1. In the 
remaining case, p > q + 1, the series diverges for i^O. The analytic 
continuation of the series q +iF q has been recently considered in [Tl][T2]. 
In this section a brackets representation of the hypergeometric series is 
obtained and then employed to produce its analytic extension. 

Theorem 5.1. The bracket representation of the hypergeometric func- 
tion is given by 

{aj + n + tj) yj (1 - b k — n + s k ) 



pF q \X 



n 

t\,...,tp 
si,—,s q 



3=1 



r(oj-) 



n 

fe=i 



r 1 - b k 



METHOD OF BRACKETS 



13 



Proof. This follows from ( 15.1 p and the representations 
(5.2) 



r(a i + n) 



r(%) 



r(%) 



r a j+ n-l e -r rfr 



E 



as well as 
(5.3) ' 



£(l-6fc-n) v fl r^ 



(oj + n + tj) 

r(%) 

1 - 6 fe - n + s k ) 



r(i-&* 



□ 



(6fe)n 

for the Pochhammer symbol. 

The bracket expression for the hypergeometric function given in The- 
orem 15.11 contains p + q brackets and p + q + 1 indices (n, tj and s&). 
This leads to a full rank system 

(5.4) a,j + n + tj =0 for 1 < j < p 

1 — bk — n + Sk = for 1 < k < q. 

of linear equations of size (p + q + 1) x (p + q) and determinant 1. For 
each choice of an index as a free variable the method of brackets yields 
a one- dimensional series for the integral. 



Series with n as a free variable. Solving (15.41) yields t* 



n 



and s 



p 



n with 1 < j < p and 1 < k < q. Rule 12.21 yields 



n 



iY-n+l-fe* 



E 



(oi), 



. , T{ aj ) a± r(l-6*) (&i)« •••(&„)« n! ' 

n=0 j=l J k=l n=0 v ' v H ' 

This is the original series representation (15. ip of the hypergeometric 
function. In particular, in the case q = p — 1, this series converges for 
\x\ < 1. 

Series with ti as a free variable. Fix an index i in the range 1 < 
i < p and solve (15 .4p to get n* = —a,i — ti, as well as £• = U — aj + a% 
for 1 < j < p, j 7^ i, and s* k = —(1 — bk) — a% — 1% for 1 < k < q. The 
method of brackets then produces the series 



E 

u 



1 I OLa 



n 



r(i - b k + at + tj) 
rn - h) 



which may be rewritten as 



(5.5) 



-x 



,-a.i 



n 



r(6 fc ) 



x q+lF p -i 



a U {1 — &fc + Oi}l<fc<g 
{1 — dj + ai}l<j<PJ^i 



-1) 



p+q-1 



X 
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Recall that the initial hypergeometric series p F q (x) converges for some 
x 7^ if and only if p < q + 1. Hence, assuming that p < q+1, observe 
that the hypergeometric series (15. 5p converges for some x if and only 
if p — q + 1. 



Series with Sj as a free variable. Proceeding as in the previous case 
and choosing i in the range 1 < i < q and then Sj as the free index, 
gives 



( , 6) [( _ ira - M -E|^ n r^ n r^ 

3 k^i 



X p F q 



{aj + 1 — bi]i<j< p 
2 — h, {1 — bk + bi}i<k<q,k^i 



X 



Summary. Assume p = q + 1 and sum up the series coming from the 
method of brackets converging in the common region \x\ > 1. Rule 
gives the analytic continuation 



9+1 

(5.7) 9¥1 F q (x) = Y,(-x) 
i=i 



n r ( Q i - 
r( aj ) 



n 



r(6 fc - ai) 



X q+lFq 



a i, {1 ~ bk + di}i<k<q 
{1 — Qj + Q-i\l<j<q+l,j^i 



for the series (15.11) . 

On the other hand, the q + 1 functions coming from choosing n or 
Si, 1 < % < q, as the free variables form linearly independent solutions 
to the hypergeometric differential equation 



(5.8) 



q+1 

n 

i=i 



n 

fc=i 



in a neighborhood of x = 0. Likewise, the g + 1 functions (15. 5ft coming 
from choosing i i; 1 < ? < q+1, as the free variables form linearly 
independent solutions to (15.81) in a neighborhood of x = oo. 



Example 5.2. For instance, if p = 2 and q = 1 then 
a, 6 



(5.9) 2 F 1 



x 



( r] -a T (b- a)T(c) (a,l-c + a 

" 1 ' T(b)T(c-a) 2 1 V 1-b + a 

h T(a-b)T(c) fb,l-c + b 

+ l J T(a)T(c-b) 2 1 V l-a + 6 
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This is entry 9.132.1 of [7]. On the other hand, the two functions 

(5.10) 2 Ft ( 

form a basis of the solutions to the second-order hypergeometric differ- 
ential equation 

( d \ ( d \ ( d \ 

(5.11) x— + a) [ x— + b) y = x— + c ) y 

\ dx ) \ dx j \ dx j 

in a neighborhood of x = 0. 



x 



X 



,1-c 



Fx 



a + 1 — c, 6 + 1 — c 



x 



6. Feynman diagram application 

In Quantum Field Theory the permanent contrast between experi- 
mental measurements and theoretical models has been possible due to 
the development of novel and powerful analytical and numerical tech- 
niques in perturbative calculations. The fundamental problem that 
arises in perturbation theory is the actual calculation of the loop inte- 
grals associated to the Feynman diagrams, whose solution is specially 
difficult since these integrals contain in general both ultraviolet (UV) 
and infrared (IR) divergences. Using the dimensional regularization 
scheme, which extends the dimensionality of space-time by adding a 
fractional piece (D = 4 — 2e), it is possible to know the behavior of 
such divergences in terms of Laurent expansions with respect to the 
dimensional regulator e when it tends to zero 

As an illustration of the use of method of brackets, the Feynman 
diagram 



(6.1) 



Po 



a i 



(i-2 



a 3 



Pi 



considered in [2] is resolved. In this diagram the propagator (or internal 
line) associated to the index a\ has mass m and the other parameters 
are P? = P 3 2 = and Pf = (Pi + P 3 ) 2 = s. The D- dimensional 
representation in Minkowski space is given by 



(6.2) 



G 



d D q 



tn°/ 2 [(Pi + q) 2 - m 2 ] ai [(P 3 - g) 2 f 2 [g 2 l a3 ' 
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In order to evaluate this integral, the Schwinger parametrization of 
(16. 2p is considered (see [S] for details). This is given by 

6.3 G = -V J H 

n? =1 r(a,) 

with H defined by 
(6.4) 

OOOOOO / •) \ / TITO 

C C c exp [x\m ) exp — s 

H= xl^xf^xf- 1 V D X , 2 2 ' d Xl dx 2 dx 3 . 

Ill (x 1 +x 2 + x 3 ) 

To apply the method of brackets the exponential terms are expanded 

as 

XlX2 „ \ _ A ( -Wni^ni S U2 X l X 2 



exp (a^m 2 ) exp ( — !— ^ s ) = V] n n2 (-l) ni m 2 

V Xi + x 2 + x 3 J ^ 

and then (16. 4p is transformed into 
(6.5) 



(xi +x 2 + x 3 ) 



n 2 ' 



oo oo oo 

ai+ni+B2-l tt2+«2-l 03-1 
X i Xn X q 



m,n 2 

Further expanding 

1 V <h r "3 T "4 T n 5 (f + n 2 + n3 + n 4 + n 5 ) 

(X! + x 2 + * 3 ) D / 2+ ^ ^ ^ 3 '" 4 ' n5 Xl X2 Xs T(f + n 2 ) 

n3,n,4,n5 y 2 ' 

and replacing into (16.51) and substituting the resulting integrals by the 
corresponding brackets yields 

n, 2m „ W2 (f + 7l 2+7l3+W4 + W 5 ) 



(6.6) H= ^^ } (-l) Bl m 



s 



x (ai + ni + n 2 + n 3 ) (a 2 + n 2 + n 4 ) (a 3 + n 5 ) . 

This bracket series is now evaluated employing Rules 12.21 and 12.31 Pos- 
sible choices for free variables are ni, n 2 , and n 4 . The series associated 
to n 2 converges for |-^-| < 1, whereas the series associated to ni,n 4 

2 

converge for < 1. The following two representations for G follow 
from here. 

Theorem 6.1. In the reqion 1-^1 < 1, 



(6.7) H = V2 ■ 2 Fi 



a,i + a 2 + 0,3 — , a 2 

D 

2 
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with rj2 defined by 



V2 = (-m 2 )' 



-a\—a2—az 



T(a 2 )T(a 3 )T (a x + a 2 + a 3 - f ) T (f - a 2 - a 3 ) 



r(f) 



Theorem 6.2. In the region |^-| < 1, 
(6.8) H = 77! • 2j P a 



ai + a 2 + % — t> 1 + a i + a 2 + a 3 — m 2 



1 + a x + a 3 - -y 



+ ??4 • 2-^1 



2 ' 



— dx — a 3 + -j 
?7i, ?7 4 defined by 

T(a 3 )r (ai + a 2 + a 3 - f ) r (f - ai - a 3 ) T (f - a 2 - a 3 ) 



m 



?7i = s 2 



-a 1 -a 2 -a 3 



-a-2 



(-m 2 y 



r (D — a\ — a 2 — a 3 ) 
r(a 2 )r(a 3 )r (ai + a 3 - § ) T (f - a 2 - a 3 ) 



These two solutions are now specialized to CL\ — Qi2 — d 3 — 1. This 
situation is specially relevant, since when an arbitrary Feynman dia- 
gram is computed, the indices associated to the propagators are nor- 
mally 1. Then, with D = 4 — 2e, the equations (16. 7p and (I6.8P take the 
form 



(6.9) H = (-m 2 ) 



< 1, as well as 



2\-l-e 



l + e,l 
2-e 



for 

(6.10) 

T(l-2e) V s 



- 2e r(e) 

2-Tl 

es 



m 
s 



for | — | < 1. Observe that these representations both have a pole at 
e = of first order (for the second representation, each of the summands 
has a pole of second order which cancel each other). 



7. Conclusions and future work 

The method of brackets provides a very effective procedure to eval- 
uate definite integrals over the interval [0, oo). The method is based 
on a heuristic list of rules on the bracket series associated to such in- 
tegrals. In particular, a variety of examples that illustrate the power 
of this method has been provided. A rigorous validation of these rules 
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as well as a systematic study of integrals from Feynman diagrams is in 
progress. 
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